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Lattice QCD at finite isospin density and/or temperature* 

J. B. Kogut^t and D. K. Sinclair, 
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^HEP Division, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, IL 60439, USA 

We simulate two-flavour lattice QCD with at a finite chemical potential fij for isospin, and finite temperature. 
At small fii , we determine the position of the crossover from hadronic matter to a quark-gluon plasma as a function 
of fij. At larger fij we observe the phase transition from the superfluid pion-condensed phase to a quark-gluon 
plasma, noting its change from second order to first order as fij is increased. We also simulate two-flavour lattice 
QCD at zero quark mass, using an action which includes an additional 4-fermion interaction, at temperatures 
close to the chiral transition on A'^j = 8 lattices. 



1. Introduction 

QCD at finite temperature and/or densities is 
relevant to the physics of the early universe, neu- 
tron stars and relativistic heavy-ion collisions — 
RHIC and the CERN heavy-ion program. 

Because of the difficulties involved with sim- 
ulating QCD at finite chemical potential fi for 
quark number with its complex fermion deter- 
minant, we are simulating Nf — 2 lattice QCD 
with a finite chemical potential fij for isospin 
(73), which has a real positive fermion determi- 
nant P]. We are currently studying the /i/ de- 
pendence of the finite temperature transition 
At small /i, the Bielefeld-Swansea collaboration 
has observed that the phase of the fermion deter- 
minant is sufficiently well behaved that the fi and 
fij dependence of the transition are identical '3] . 
Our predictions are in good agreement with those 
of de Forcrand and Philipsen j^. 

At large values of /i/ (/x/ > m.^), the low tem- 
perature phase is characterized by a pion conden- 
sate. The finite temperature transition is now a 
true phase transition, which appears to be second 
order for lower fj,i values and first order at higher 
HI values. 

We also simulate lattice QCD with an irrele- 
vant chiral 4-fermion term which allows us to sim- 
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ulate at zero quark mass, giving us direct access 
to the critical exponents at the finite temperature 
transition 0. We are using 16^ x 8 and 24^^ x 8 
lattices for these simulations. 

In section 2 we present preliminary results of 
our finite /i/ and temperature simulations. Sec- 
tion 3 gives some preliminary graphs from our 
Nt — S finite-temperature simulations using our 
modified action. Finally we present our conclu- 
sions and indicate further avenues of research in 
section 4. 

2. QCD at finite fij and temperature 

The quark part of our lattice action is 
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(1) 



For simulations at /i/ < m^r, we set the symmetry 
breaking parameter A = 0. 

We perform simulations on an 8^ x4 lattice with 
m — 0.05, A = and < ^/ < 0.55, for a set of 
(3s covering the crossover region for each fii. We 
measure the chiral condensate, the plaquette, the 
Wilson Line and the isospin density and their cor- 
responding susceptibilities for each set of param- 
eters. The position of the crossover was obtained 
as the peak of the susceptibilities. This was de- 
termined using Ferrenberg-Swendsen reweighting 
techniques 0. Figure shows the chiral suscep- 
tibilities obtained from such reweightings for 3 of 
the 7 /i/ values we use. The multiple values for 
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each Hi are the resuhs of using distributions from 
several /3s close to the peak. 

SU(3) N,=2 m=0.05 X=0 8^x4 lattice 
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Figure 2. Position of the crossover as functions of 



Figure 1. Ferrenberg-Swendsen reweighted chiral 
susceptibilities for fii = 0.0, 0.4, 0.55. 



The crossover f3 values, (3c, obtained from the 
peaks in the susceptibilities for each of the 4 
observables, are plotted in figure 13 as functions 
of /ij, since the leading term is expected to be 
quadratic in fij. We note that the 4 predictions 
for each fii appear consistent within errors, indi- 
cating that this is a reasonable definition of the 
position of the crossover. The straight line 



f3c = 5.322 



(2) 



in this figure is only meant as a rough guide. 
Using the relation /3c (m) — 0c{f^i = 2^), which 
should hold for small /i, /i/ our results are consis- 
tent with those of de Forcrand and Philipsen 0]. 

We are repeating these simulations for m = 0.1 
and m = 0.2 where the critical fxjs, fic — "t-tt 
are larger. In no case does there appear to be 
a critical end point, beyond which the transition 
becomes first order, for fij < m.^^. 



We are extending our m — 0.05 simulations to 
fii > m^r. While the thermal phase transition 
where the pion condensate evaporates appears to 
be second order for fij = 0.6, simulations on 16"^ x 
4 lattices show indications of the metastability 
expected for a first order transition, at /z/ = 0.8. 



3. "xQCD" at finite temperature 

We simulate lattice QCD with 2 flavours of 
massless staggered quarks and an irrelevant chi- 
ral 4-fermion interaction which allows us to run at 
zero quark mass, at finite temperature on 16'^ x 8 
and 24^ x 8 lattices. Using Nt = 8 should free us 
from the lattice artifacts which were present for 
A^t = 4,6 0. 

Having massless quarks should enable us to ex- 
tract the critical exponent /?„ which describes 
the vanishing of the chiral condensate as the chi- 
ral transition is approached from below, and the 
critical (3 (/3c)- Running at /3c with small quark 
masses will give us the critical exponent S. 

Preliminary results for the chiral condensate 
and Wilson Line are given in figure |31 These 
clearly indicate that the transition occurs in the 
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range 5.530 < Pc < 5.545 and appears sufficiently 
smooth to be second order. 
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Figure 3. Chiral condensate and Wilson line of 
"xQCD" as functions of /3 on A^j = 8 lattices. 



4. Conclusions 

We have determined the fj,j dependence of the 
finite temperature transition for 2-flavour lattice 
QCD. For small ^, the fluctuations of the phase 
of the fermion determinant are small enough that 
the and /z/ dependence of this transition are 
the same [Sj- The /Lt/ dependence we measure 
predicts a n dependence consistent with that ob- 
tained from simulations with imaginary /j, 0]. We 
will check if it is also consistent with series ex- 
pansions around ^ = Our simulations at 
3 different quark masses indicate that the tran- 
sition is a crossover for /i/ < m.^-. At higher fij, 
where the the pion condensate evaporates at the 
transition, we see evidence for a change to first 
order behaviour, but it is unclear if this is related 
to the critical endpoint expected at finite fj,. 

Others have estimated the position of the crit- 
ical endpoint at finite /i for the more physical 
2-|-l-flavour QCD [3, the closely related 3-flavour 
QCD Hi and 4-flavour QCD We are ex- 



tending our simulations to 3-flavours, where one 
can tunc the quark mass to make the critical end- 
point as close to /i = as desired. If the endpoint 
is at small enough /x, we expect a critical endpoint 
at the corresponding /z/, giving another estimate 
for the position of this endpoint. 

We are using the "xQCD" action to determine 
the critical exponents for the chiral transition of 
finite temperature QCD, on A^t = 8 lattices. 

We thank F. Karsch, P. de Forcrand, 

0. Philipsen and D. Toublan for helpful discus- 
sions. These simulations were performed on the 
IBM SP at NERSC and the Jazz cluster at Ar- 
gonnc. 
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